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Abstract 

We consider the motion of test particles in the spacetime of a black hole in f{R) gravity. The 
complete set of analytic solutions of the geodesic equation in the spacetime of this black hole 
are presented. The geodesic equations are solved in terms of Weierstrass elliptic functions and 
derivatives of Kleinian sigma functions. The different types of the resulting orbits are characterized 
in terms of the conserved energy and angular momentum as well as the cosmological constant (A) 
and the real constant (/3). 
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1. INTRODUCTION 


Recent observations of the Supernova Type la (SNIa) and Cosmic Microwave Background 


Q. 


(CMB) radiation indicate that the expansion of our Universe is accelerating [ll. This ac¬ 
celerating expansion is one of the most important puzzles of contemporary physics. While 
a nonzero vacuum energy could cause the Universe to accelerate, one is left wondering why 
the vacuum energy would be so small { 2 , 3|. As a consequence, all observations related to 
gravity should be described within a framework including the cosmological constant. How¬ 
ever, due to the smallness of the cosmological constant it seems unlikely that it will have an 
observable effect on smaller, that is, on solar system scales. In fact, it has been shown that, 
based on the Schwarzschild-de Sitter spacetime, the cosmological constant plays no role in 
solar system observations 4| and neither in strong helds j^. 

Another question concerning the cosmological constant is the coincidence problem, i.e., 
the fact that the energy density of matter and of the vacuum (cosmological constant) are 
on the same order of magnitude at the present time. A model of varying dark energy can 
partially solve this problem. Here the density of dark energy traces the density of matter 
from the early universe to the present time. Modified gravity can also provide an effective 
time varying equation of state. In such models the Einstein-Hilbert action is replaced with 
a generic form of f{R) gravity But modifying the action not only affects the dynamics 
of the Universe, it can also alter the dynamics at the galactic or solar system scales. 

Although the majority of gravitational effects can be discussed using approximations and 
numerics, a systematic study of all effects can only be achieved by using analytical methods. 
The motion of test particles (both massive and massless) provides the only experimentally 
feasible way to study the gravitational helds of objects such as black holes. Predictions 
about observable effects (such as light dehection, gravitational time-delay, perihelion shift 
and the Lense-Thirring effect) can be made and compared with observations. For this 
purpose the Weierstrassian elliptic functions are most useful because they lead to simple 
expressions. The resulting structure of the equations of motion is essentially the same as 
in the Schwarzschild spacetime, where they can be solved analytically in terms of elliptic 


Q. 


functions, as hrst demonstrated by Hagihara in 1931 

This analytical method was further advanced and applied to the hyperelliptical case, 
where the analytical solution of the equations of motion in the 4-dimensional Schwarzschild- 
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de Sitter and Kerr-de Sitter spacetime [8[ was presented. It was applied as well as in the 
higher-dimensional Schwarzschild, Schwarzschild-(anti)de Sitter, Reissner-Nordstrom and 
Reissner-Nordstrom -(anti)-de Sitter spactime and in the higher-dimensional Myers-Perry 
spacetime ^ 11|. Also the general solutions to the geodesic equation in Schwarzschild and 
Kerr black hole pierced by a cosmic string have been discussed extensively [3, Q- Recently, 
the geodesic equations were solved analytically in special cases in the singly spinning black 


ring spacetime [l^, as well as in the (charged 


and in the (rotating) black string spacetime 16| . 


doubly spinning black ring spacetime 


id 


In this paper we discussed the motion of test particles in the spacetime of the black hole 
derived in modihed gravity . We present the results here in terms of Weierstrass elliptic 
functions and derivatives of Kleinian sigma functions. Our paper is organized as follows: in 
Section ([2]) we give a brief review of the held equations in f{R) modihed gravity. In Section 
([3]) we present the geodesic equations describing test particle motion in the spacetime of this 
black hole and we give our analytical solutions. We conclude in Section (j4]). 


2. FIELD EQUATIONS IN f{R) MODIFIED GRAVITY 


In this section we give a brief review of the held equations in f{R) gravity. A generic 
form of the action depending on the Ricci scalar can be written as follows: 




F{R)Rf,^ - ^f{R)g^^ - - gf,uO)F{R) = 


Varying the action with respect to the metric results in the held equations as: 

1 
2 

where F{R) = ^ and □ = V„V“. 

The metric of the spherically symmetric spacetime we are considering is given by [1 

ds^ = B{r)dt^ — A{r)dr'^ — r‘^{d6‘^ + sin^ Od^), 

, where A{r) = . The model employed for f{R) gravity is given by 


( 1 ) 


( 2 ) 


f{R) = R +A + 


R +A 


In 


R +A 


( 3 ) 


( 4 ) 


R/Rq -|- 2/q; R^ 

where A is the cosmological constant, which has a value of |A| < j^, Rc is a 

constant of integration and Rq = where a and d are free parameters of the action. 
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The metric solution up to the hrst order in the free parameters of the action is obtained as 
i?(r) = 1 — ^ + /3r — |Ar^, where (3 = a/d > Q is a. real constant 17|. In this model d is 
a scale factor around 10 kpc, and a ~ 10“® is a dimensionless parameter which provides a 
flat rotation curve for stars in a typical spiral galaxy. 

In the case of strong curvatures where A and R/Ro 3> 2/a, the action reduces to 

R 


f{R) = R + Roin 


Rr' 


( 5 ) 


which is the desired action for the solar system, typical black holes as well as spiral galaxies. 
On the other hand on cosmological scales, where i? ~ i?o — and when a <C 1, the action 
reduces to f{R) = R +A which satishes the late time accelerating expansion of the Universe 


izl 


Therefore the small free parameter, which is denoted by (3 in this paper, satishes both 
the small and the large scale physical evidence for the two ranges of the scalar curvature, 
S> A and i? ~ A. 

3. THE GEODESIC EQUATION 


We consider the geodesic equation 

d3^x^ dxP dx'^ 

-L -= 0 

ds ds 

where ds^ = g^ydx^dx'' is the proper time along the geodesics and 

Tpa = ^g^^idpOau + d^gpy - dygp^) 
is the Christoffel symbol, in a spacetime given by the metric 


( 6 ) 


( 7 ) 


ds^ = {1 -^ +/3r —-Ar^)(it^ — (1 — +/dr —-Ar^) ^dr"^ — r‘^{d9‘^+ siv? 9dp^), (8) 

r 3 r 3 

which describes the spherically symmetric vacuum solution of Eq.([T]). For a general dis¬ 
cussion of this metric, see e.g. 17|. The Lagrangian £ for a point particle in the spacetime 
Eq.Q reads : 

^ 1 dxP dx'^ 1 Ir^r^/ ..^dt.^ , T,dr.^ ^,,d9,^ . o^,dip.n,, 

^ ) > 1 - 


where for massive particles e = 1 and for light e = 0 
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Because of the spherical symmetry we can restrict our considerations to the equatorial 
plane. Furthermore, due to the conserved energy and angular momentum 


^ dt 2m ^ 1 ^ 2 \ 

ds^ 

the geodesic equation reduces to an ordinary differential equation 

,dr 


ds' 


r 




( 10 ) 

( 11 ) 

( 12 ) 


Together with energy and angular momentum conservation we obtain the corresponding 
equations for r as a function of 0 and as a function of t 


,dr.r, , 2 m ^ 1 . 9 ,, , 

=J 2 (^ -{1- —+ (3r--Ar){e+ —))=: R{r), 


(13) 


'—R = 
'dt^ 


2m ^ 1 . 9 x 9 / 9 / 2 m 

-h /Sr — -AP) (E^ — (1-h (3r 

T o T 


iArX(£ + ^)). (14) 


Eqs. (IT^ - flIT|) give a complete description of the dynamics. Eq.([T2]) suggests the intro¬ 
duction of an effective potential 

S^7T7 1 T 

^e// = (l- —+ /?r--Ar^)(e + -). (15) 

The shape of an orbit depends on the energy E and the angular momentum L of the test 
particle or light ray under consideration as well as the cosmological constant A and the real 
constant fd. (The mass can be absorbed through a rescaling of the radial coordinate.) Since 
r should be real and positive, the physically acceptable regions are given by those values of 
r, for which E"^ > 14//, owing to the square on the left hand side of Eq. (ll2p . 

The following different types of orbits can be identihed in the spacetimes described by 
the metric Eq.([H]) 

1. Flyby orbits (FO): r starts from cx), then approaches a periapsis r = Vp and goes back 
to CX). 

2. Bound orbits (BO): r oscillates between two boundary values Vp < r < Va with 
0 < Tp < Tq < CX) . 
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3. Terminating bound orbits (TBO): r starts in (0,ra] for 0 < < oo and falls into the 

singularity at r = 0. 

4. Terminating escape orbits (TEO): r comes from cxo and falls into the singularity at 
r = 0. 

The four regular types of geodesic motion correspond to different arrangements of the 
real and positive zeros of R{r) dehning the borders of R{r) > 0 or, equivalently, E'^ > Vg/f . 

All other types of orbits are exceptional and treated separately. They are connected with 
the appearance of multiple zeros in R{r) or with parameter values which reduce the degree 
of R{r). In both cases the structure of the differential Eq. flldp is considerably simplihed. 
These orbits are radial geodesics with L = 0 (i.e., ^ = 0), circular orbits with constant 
r, orbits asymptotically approaching circular orbits, and in the case of /9 = 0 and A = 0 
parabolic orbits with E‘^ = 1 . 

For the analysis of the dependence of the possible types of orbits on the parameters of the 
spacetime and the test particle or light ray it is convenient to use dimensionless quantities. 
Thus, we introduce 


T ~ 1 ~ 

r = —, £ = —, A = -Am^ (5 = /3m. (16) 

m 6 

Eq. ffT3|) can be written as 

(I'T ~ ~ ~ ~ 

( —= eACr^ - e/3Cr^ + {{E^ - e)C + A)r^ + (2e£ - /3)f^ -r^ + 2r = R{r). (17) 

dip 

Eq. flTT)) implies that R{r) > 0 is a necessary condition for the existence of a geodesic and, 
thus, that the real and positive zeros of R{r) are extremal values of the geodesic motion. In 
the limit /3 = 0,A = 0, Eq. flTTl) reduces to the well known Schwarzschild limit, and in the 
limit /3 = 0 to the Schwarzschild-de Sitter limit. As r = 0 is a zero of R{r) for all values of 
the parameters, it is neglected in the following analysis and 

R*{r) = eACf^ — + ((E^ — e)C + A)r^ + (2e£ — /3)r^ — r + 2 (18) 

is considered instead. 

If for a given set of parameters A,/3,e,E‘^, and C the polynomial R*{f) has n positive 
and real zeros, then for varying and L this number can only change if two zeros merge 
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to one. Solving R*{r) = 0, = 0 for E'^ and C, for e = 1 yields 


dr 

:t3 _ /5^2 


= 


and for e = 0 


C = 


2(Ar^ — — r + 2)^ 

f(/3f^ + 2f — 6) 

T , 4/3^ 


C = 


-{i3r‘^ + 2r — 6) 
f2(2Af3 — /3f2 — 2) 


^2 


-A). 


( 19 ) 


( 20 ) 


(-l + \/l + 6 /?)i> (- 1 +\/l + 6 /l )2 

In the limit /3 —)■ 0 Eqs. flT^ . fl20l) rednce to the Schwarzschild-de Sitter expressions = 


(Af®—r+2)^ ^ _ —(f—3) 


and for e = 0 , £ = — A) and in the limit /3 —)• 0 and A —)■ 0 


f(r-3) ’f2(Ar3-l) '' ^2 V 27 

reduce to the Schwarzschild expressions E^ = 


and for e = 0 , £ = 27 ^ 


In 


Figs. ffT]|2]l3] and I4|) the results of this analysis are shown for both test particles (e = 1 ) and 
light (e = 0) for a small positive cosmological constant A and /3 . 
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(a) Schwarzschild space-time with A=0, |3=0. 



(b) Schwarzschild space-time with A=0, (3=10“^ 


FIG. 1: Regions of different types of geodesic motion for test particles (e = 1). The numbers of 
positive real zeros in these regions are: (I)=2, (II)=0, (III)=3, (IV)=1. 
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(c) Schwarzschild space-time with A=0, (3=10 ^ 



£■2 


(d) Schwarzschild-de Sitter space-time with A= "^ (10 ^ 

FIG. 2: Regions of different types of geodesic motion for test particles (e = 1). The numbers 
of positive real zeros in these regions are: (I)=2, (II)=0, (III)=3 for (Schwarzschild) and 4 for 
(Schwarzschild-de Sitter), (IV)=1. 
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0.12 




(1) spacetime of a black hole in f(R) gravity with A=-^ (10 ^), P=10 ^ 


FIG. 3: Regions of different types of geodesic motion for test particles (e = 1). The numbers of 
positive real zeros in these regions are: (I)=2, (II)=0, (III)=3. For the effective potentials see 
Fig. dSD 
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(g) Null, £ =0. 

FIG. 4: Regions of different types of geodesic motion for light (e = 0). The numbers of positive 
real zeros in these regions are: (I)=2, (II)=0. For the effective potentials see Fig. (I5|) 
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Three different regions of geodesic motion can be identified: 

1. R*{r) has 2 positive real zeros ri < r 2 with R{r) > 0 for 0 < r < ri and r 2 < r. 
Possible orbit types: flyby and terminating bound orbits. 

2. R*{r) has 0 positive real zeros and R{r) > 0 for positive r. Possible orbit types: 
terminating escape orbits. 

3. R*{r) has 4 positive real zeros r* < rj+i with R{r) > 0 for 0 < r < ri, r 2 < r < rs and 
r 4 < r. Possible orbit types: flyby, bound, and terminating bound orbits. 

TABLE I: Types of orbits of light and particles in the spacetime of a black hole derived from the 
f{R) modihed gravity. The second column gives the number of positive zeros of the polynomial 
R*(r). In the third column, the thick lines represent the range of orbits, and the turning points 
are shown by thick dots. The small vertical line denotes f = 0. 


Region 

Pos. Zeros 

Range of r 

Taype of Orbits 

I 

2 

-• --> 

(TBO) , (FO) 

II 

0 

-> 

(TEO) 

III 

4 

i— • — • — • — • — > 

(TBO) , (BO) , 

__ 


For light rays only regions (I) and (II) appear. Sample effective potentials for each of 
the regions (I), (II), and (III) are shown in Fig. ([5]). A summary of possible orbit types 
can be found in Tab. ([T]). For light as well as test particles exceptional orbits appear at 
the boundaries of the regions (I) to (HI), corresponding to multiple zeros of R*{r). In the 
case of multiple zeros the boundary is described by and C given by Eq. flT^ or fl20|) . 
respectively, and the corresponding test particle or light ray moves on a circular orbit, which 
may be stable or unstable. If R*{r) has a maximum, then the circular orbit is stable and 
else unstable. For test particles, the substitution of Eq. flT^ in yields 


d^R*{f) _ -2(3A/3f^ + (8A - - 30Ah=^ + 12 ^^^ _ 2 y + 12 ) 

r‘^(2Ar^ — — 2 ) 

In the limit = 0 Eq. fl2lll reduces to the Schwarzschild-de Sitter expression 
-2{4Ar^ ^-i^5Ar ^-r+6) limit /3 = 0 and A = 0 to the Schwarzschild expression 

-2(f-6) 

f.2 


( 21 ) 


d?R*{f) 

df'^ 

df^ 
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The boundary of region (III) is defined by two corners, which correspond to the triple 
zeros of R*{r). This means that the region will vanish if triple zeros are not possible. With 
an analysis similar to Eq. fll9p it can be shown that triple zeros appear at 

^2 8(/3f^ + 2r — 6)^ 

~ h(3/3f2 + 8h-30)2’ 

^ _ (3/3^2 + 8h - 30) 

P0P + 6) 

^ + 3/3h3 - 12/3f^ + 2h - 12) 

P{3pp + 8r-30) 


In the limit (3 = 0 Eqs. fl2^ reduce to the Schwarzschild-de Sitter expressions = 

16(f-3)^ ^ _ 4f-15 t _ (r-6) 

f(4f—15)^ ’ 3r^ ’ f^(4'r—15) * 

Let us now compare the geodesic motion of test particles (e = 1) with the case (3 = 0. 
From Figs. flT]|20 and 0]) it is obvious that some regions for small positive (3 are a little 
shifted compared to /? = 0. But for the larger positive /3, all regions for (3 > 0 are shifted 
compared to /3 = 0. Thus, every pair £), which was located near a boundary for (3 = 0, 
may switch to another region. Regions (IV) and (III) absorbed the whole parameter space, 
which was located in regions (II) and (I), respectively, for /3 = 0, see Fig. (1T1I21 a,c) and 
for cases with the positive cosmological constant, regions (I) absorbed the whole parameter 
space, which was located in regions (II) for (3 = 0. For (3 = 0, there are also terminating 
escape orbits which, coming from inhnity, are reflected at the potential barrier induced by 
the larger positive (3. The types of orbits for null geodesics, i.e., for light rays, are not 
essentially influenced by a nonvanishing (3. For every region, examples of timelike and null 
geodesics can be found in Figs. fl6H7ll8]IMT0]llllll2l) . 


3.1. Analytical solution of geodesic equations 

In this section we present the analytical solution of the equations of motion. We introduce 
a new variable u = I, and obtain from Eq. flTT)) : 

(^)2 = 2u^ -u^ + (2eC - P)u + {{E^ - e)C + A) - e^C- + ekC^. (23) 
(icp u 
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3.1.1. Null geodesics 


For e = 0 Eq. (1^ is of elliptic type 

,dU,r 


'dip' 


2u^ — + {E^C + A) = P3 (m) = 


y ^ dill 
i=0 


(24) 


A further substitution u = ^(4|/ — ^) = 2?/ + ^ transforms pziu) into the Weierstrass form 


so that Eq. fl2^ turns into: 


(^)2 = 4^3 _ 

dip 


(25) 


where 


, aiOaUg aoaj a| /3 {E^C + A) 1 

Qo — — - — —, Q"^ — - — - — - — — - ~h - (■^o) 

^ 12 4 12 2 ’ ^ 48 16 216 24 4 216 ^ ^ 


al aiOg 1 (3 


are the Weierstrass invariants. The differential equation fl2^ is of elliptic type and is solved 
by the Weierstrass function 0 ,H,Q- 


y{l^) = p{^ - Pin] 92,93), 


(27) 


where ipin = ipo + 


dy 


yo 


\/ ‘iy^-92y-. 


93 


’^ith j/o = || + f| = l(T-i), depending only on the 


initial values ipo and tq. Then the solution of flTTl) acquires the form 

03 1 


r{p) = 


3 ‘^p{p - Pin] 92, 93) + Q 


Ap{p - Pin] 92, 93 ) - f 
The corresponding light trajectories have been exhaustively discussed in 0. 


(28) 


3.1.2. Timelike geodesics 


For e = 1 Eq. fl2^ should be rewritten as 


du , 


{u——Y = 2u^ — u‘^ + {2C — /3)u^ + ((E^ — 1)£ + A)u^ — ^Cu + AC = pAu) 

dip 


udu 


(29) 


i=0 


A separation of variables in fl29|) yields 


p-Po = 


'no 


(30) 
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where uq = u^tpo). In solving integral (l30|) there are two major issues which have to be 
addressed. First, the integrand is not well dehned in the complex plane because of the two 


branc 

path 


les of the square root. Second, the solution should not depend on the integration 
If 7 denotes some closed integration path and 

udu 


CO = 


this means that 


Vp5(m)’ 

/■“ udu 


if- ipo-u = 


should be valid, too. Hence, the solution u{ip) of our problem has to fulhll 

u{ip) = u{ip — cj) 


(31) 


(32) 


(33) 


for every ca 7 ^ 0 obtained from Eq. fl3Tl) . These two issues can be solved if we consider Eq. 
(130]) to be defined on the Riemann surface = p^i^x) of genus g = 2 and introduce a basis 
of canonical holomorphic and meromorphic differentials dzi and dvi, respectively. 


dzi = 


dx 




dzo = 


xdx 


a^x + 2a4X^ + Sa^x^ 

(IT I = - ^ -ax, 


dro = 


x'^dx 


4\/p^ ’ ^ 4^p5(a;)’ 

and real 2cjij , 2pij and imaginary 2cj[j , 2p[j period matrices 


Now we can write the analytic solution of equation 


(34) 

(35) 



f dzi, 

24 = 

f dzi, 

( 36 ) 


Jaj 




(h dvi, 

H,= 

® dvi. 

( 37 ) 


as 


u{p) = - —((p^), 

<72 

where Uj is the i-th derivative of the Kleinian sigma function in two variables 

a{z) = Ce^"’^^e[Kj\{2u-^z-T), 

which is given by the Riemann 0-function with characteristic RToo, 

e{z-,T) = V 


( 38 ) 


( 39 ) 


( 40 ) 
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A number of parameters enter here: the symmetric Riemann matrix r = the 

period-matrix (2a;, 2a;'), the period-matrix of the second kind {2ri,2ri'), the matrix k = 
t]{2u)~^, and the vector of Riemann constants with base point at inhnity 2Koo = (0, 1 )* -|- 
(1, l)'r. The constant C can be given explicitly, see e.g. [2^, but does not matter here. 
In Eq. 0381) the argument is an element of the one-dimensional sigma divisor: = 

{f{ip - p>in), if - ifinY where ipin = (^0 + /“ Uq = ^ depends only on the initial 

values, and the function / is given by the condition = 0. For more details on the 


construction of such solutions see e.g.[ld, l20|. With Eq. (l38|) the solution for r is given by 


r = - — (41) 

In this equation, a\ and a 2 depend on the set of parameters ipa, 0 J,r],T and also on P 5 {u) 
according to Eqs. (1301— I 37p . which contains the /9-dependence of the modihed gravity. The 
solution of r is the analytic solution of the equation of motion of a test particle in the 
spacetime of a black hole derived in f{R) modihed gravity. This solution is valid in all 
regions of this spacetime. 


3.1.3. The orbits 

In Figs. fl6]l7ll8]IMT0]|lllll2p some of the possible orbits are plotted. For all orbits in each 
region of the hgure the parameters E and C are the same. The value of the cosmological 
constant is chosen as A = 1(10“^) in all plots. For massive test particles: In Fig. ([6]) the 
parameters = 1.1 and C = 0.025 belong to the (I) region in Fig. (|3])(e) and Fig. (I5])(a). For 
a positive cosmological constant this corresponds to a hyby orbit and a bound terminating 
orbit, ending in the singularity. The corresponding orbits are shown in Figs. l|6|)(a),(|6|)(b). 
The next parameter choice is = 2 and C = 0.11. These parameters he in the (II) region 
of Fig. (I3])(e) and Fig. (I5])(b) corresponding to a terminating escape orbit ending in the 
singularity. Fig. ((Tj). Our third choice of parameters is E‘^ = 0.93 and C = 0.072. These 
parameters are in the (III) region in Fig. (I3])(e) and Fig. (I5])(c) with four zeros, indicating a 
terminating bound orbit, a bound orbit and a flyby orbit, see Figs. (IZ])(a),(l8])(b) and (IH])(c). 

The dependence of the perihelion shift of a bound orbit of a massive test particle is shown 
in Fig. ([8])(b). The periastron advance for such a bound orbit is given by the difference of 
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icity of the solution r((^). For more details 


22| . The orbit in Fig. ©(c) is a reflection 


the 27r-periodicity of the angle ip and the perioc 
on the construction of such solution see e.g. js, 
at the A barrier. 

For light rays: In Figs. fl^llTUD some of the possible orbits are plotted. The parameters 
= 0.5 and C = 0.025 belong to the (I) region in Fig. ©(b) and Fig. ©(a), corresponding 
to a flyby orbit and a bound terminating orbit ending in the singularity. The corresponding 
orbits are shown in Figs. ([9])(a),([9])(b). The next parameter choice is = 2 and C = 0.11. 
These parameters lie in the (II) area of Fig. (ll])(b) and Fig. ([5])(b), corresponding to a 
terminating escape orbit ending in the singularity, Fig. ffHJ]) . The deflection of light (massless 
test particle, i.e. e = 0) is shown in Fig. (in])(b). For more details on the deflection of light 
solution see e.g. j^, 22|. Some of the types of orbits for timelike geodesics which are 
essentially influenced by a nonvanishing [3 are presented in Figs. fllllll2p . 
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(a) Region I 
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(b) Region II 
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(c) Region III 


FIG. 5: Effective potentials for different regions of geodesic motion. The horizontal grey dashed 
line represents the squared energy parameter 
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FIG. 6: Timelike geodesics for the different regions of orbit types (A = |(10 ^),/3 = 10 ^ 
always indicate the Schwarzschild radius. 


). Circles 
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FIG. 7: Timelike geodesics for the different regions of orbit types (A = |(10 ^),/3 = 10 ^). Circles 
always indicate the Schwarzschild radius. 
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FIG. 8: Timelike geodesics for the different regions of orbit types (A = |(10 ^),/3 = 10 ^). Circles 
always indicate the Schwarzschild radius. 
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FIG. 9: Null geodesics for the different regions of orbit types (A = ^(10 ^),/3 = 10 ^). 
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FIG. 10: Null geodesics for the different regions of orbit types (A = ^(10 ®),/3 = 10 ^). 
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FIG. 11: Comparison between orbits in the cases of nonvanishing /3 and /3 = 0. 
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FIG. 12: Comparison between orbits in the cases of nonvanishing /3 and /3 = 0. 
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4. CONCLUSIONS 


In this paper we considered the motion of massive and massless test particles in a black 
hole metric of f{R) gravity, presented in 17|. After reviewing the spacetime and the corre¬ 
sponding eqnations of motion we classihed the complete set of orbit types for massive and 
massless test particles moving on geodesics. In particular, we analyzed conditions for circu¬ 
lar orbits. The geodesic equations can be solved in terms of Weierstrass elliptic functions 
and derivatives of Kleinian sigma functions. 

We also considered all possible types of orbits. Using effective potential techniques and 
parametric diagrams, the possible types of orbits were derived. For null geodesics FO, TBO 
and TEO are possible, while for timelike geodesics BO, FO, TBO and TEO are possible. 
We also performed a comparison between orbits in the cases of nonvanishing /3 and (3 = 0. 

The results obtained in this paper present a useful tool to calculate the exact orbits 
and their properties, including observables like the periastron shift of bound orbits, the 
light deflection of flyby orbits, the deflection angle and the Tense-Thirring effect. For the 
calculation of the observables analogous formulas to those given in js, 21, 22| can be used. 
It would be interesting to extend this work to a charged and rotating version of the black 
hole spacetime, which has still to be derived. 
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